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The charge asymmetry of W bosons produced in pp collisions at i/s = 1.96 TeV is sensitive to the ratio of d 
and u quark distributions in the range of x > 0.002 at Q 2 w My^, We propose an analysis technique to directly 
measure W production charge asymmetry from W — > ev events at the Tevatron and show the feasibility for 
this method using Monte Carlo simulations. 



I. INTRODUCTION 

The differential cross section for W boson production in pp 
as a function of W rapidity is 

da^~ 2tt Gp ^-^\ . .2 r , , _, x _, 

~ A = "7 — 7= >. [q{x p )q(xp) + q(x p )q(xp)\ , 

dyw >> V2^ 

(1) 

where x p (xp) is the fraction of the proton (anti-proton) mo- 
mentum carried by the struck quark, q and q are the quark and 
anti-quark parton distribution functions, and yw is the rapid- 
ity of the W boson. The x values of the quark in the proton 
and antiquark in the antiproton are related to the rapidity, y, of 
the W boson via the equation x Pt p = Mw /\/se ±yw as shown 
in FigQ] Here y/s is the center of mass energy and Mw is the 
mass of the W boson. 

Since W + (W~ ) bosons are produced in pp collisions pri- 
marily by the annihilation of u(d) quarks in the proton and 
d(u) quarks in the anti-proton, and since u(x p ) = u(xp) and 
d(x p ) = d(xp) by CPT symmetry, the differential cross sec- 
tions for are approximately 
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FIG. 1 : The x values of the quark for W production at the Tevatron. 
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FIG. 2: (a) The W boson and lepton rapidity distributions in pp 
collisions, (b) The charge asymmetry for W production and the 
decay lepton. 



Since the u quark tends to carry a larger fraction of the pro- 
ton's momentum than the d quark on average, the 
is boosted in the proton (anti-proton) direction as shown in 
Fig. |2(a)| The W production charge asymmetry, A(yw), in 
the leading-order parton model is therefore 



A(y w ) 



da + 1 dyw — da jdy 



w 



da + 1 dyw + da ~ /dyw 
u{xp)d(xp) — d{xp)u{xp) 
u(x p )d(xp) + d{xp)u{xp) 

Rdu{xp*) Rdui.Xp) 



(4) 



Rdu(xp) + Rdu(x p ) ' 

where we use Eq.|2]and Eq.[3]and introduce the ratio Rd u = 
^2y. As we see in Eq. |4j there is a direct correlation between 
the W production charge asymmetry and the d/u ratio. A 
precise measurement of the W production charge asymmetry 
serves as a valuable constraint on the u and d quark momen- 
tum distributions [Q]]. 

W production at hadron colliders is identified through the 
process p + p — * W ± ; — > + v. Since the W de- 
cay involves a neutrino whose longitudinal momentum is ex- 
perimentally undetermined, the charge asymmetry previously 
has been constrained by the measured charge asymmetry of 
the decay leptons and as a function of the lepton pseudo- 
rapidity IJ, yl Bfl . However, as shown in Fig. |2(b)| there is a 
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"turn-over" in the lepton charge asymmetry due to a convolu- 
tion of the W production charge asymmetry and the W V — A 
decay. This convolution means leptons from a single pseudo- 
rapidity come from a range of W rapdities and thus a range 
of parton x values. Thus, the measured lepton asymmetry is 
more complicated to interpret in terms of quark distributions, 
and we expect the direct measurement of the asymmetry of 
the rapidity distribution to be a more sensitive probe of 
the differences between u and d quarks. 

In this paper, we propose a new analysis technique which 
resolves the kinematic ambiguity of the longitudinal momen- 
tum of the neutrino to directly reconstruct the rapidity. 
We describe the details of our new analysis technique and 
outline the sources of systematic uncertainty of this measure- 
ment. Our studies are performed in the — + e^v channel 
produced in pp collisions at the Tevatron. We use a realistic 
Monte Carlo simulation (MC@NLO) and include the effects 
of higher-order QCD corrections lUt] . 

II. ANALYSIS TECHNIQUE 

The W decay to leptons, in our case — ► e^v, involves 
a neutrino whose longitudinal momentum cannot be experi- 
mentally determined. However, we can determine the longitu- 
dinal momentum by constraining the W mass in Eq. [5] which 
results in a two-fold ambiguity. This ambiguity can be partly 
resolved on a statistical basis from the known V — A (vector- 
axial vector) decay distribution using the center-of-mass de- 
cay angle between the electron and the proton, 6*, and from 
the W + and W~ production cross-sections as a function of 
W rapidity, da^ /dyw- The W mass constraint is 

M% = {Ei + E v f - {Pi + P„f, (5) 

where the W mass, Mw, is contrained to its experimentally 
measured value J6l 0]. Events which cannot satisfy the W 
mass constraint (and which get imaginary values of the neu- 
trino z-momentum) are due to a mis-reconstruction of the neu- 
trino (missing) transverse energy, i£ T JH]. Therefore, in such 
cases, we re-scale the |5 T to the value which makes the imag- 
inary part to be zero. This new $ T is then used to correct the 
yw for the event. 

The leading order W boson production mechanism in pp 
collisions results in the W boson being polarized in the p di- 
rection by means of the V — A structure of the weak interac- 
tion. The V — A structure means that the weak current cou- 
ples only to left-handed u and d quarks (or to right-handed 
u and d quarks). For ultra-relativistic quarks, where helicity 
and chirality are approximately equivalent, this results in full 
polarization of the produced W bosons in the direction of the 
beam. The W leptonic decay process also couples only to 
left-handed e~ and right-handed v (or right-handed e + and 
left-handed v). The conservation of angular momentum fa- 
vors a decay with the final state lepton (neutrino or electron) 
at a small angle with respect to the initial state quark direction 
(and a similar small angle between the initial state anti-quark 
and final anti-lepton). The systematic shift in lepton pseudo- 
rapidity with respect to yw depending on the charge of the 
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FIG. 3: (a) The positively charged W boson and lepton rapidity dis- 
tribution, (b) The negatively charged W boson and lepton rapidity 
distribution. 



final state lepton is illustrated in Fig. |3(a)| and |3(b)| which 
shows the lepton pseudo-rapidity vs. W rapidity for the dif- 
ferent charges. This effect also explains the discrepancy at 
high rapidity between the lepton charge asymmetry and the 
W charge asymmetry as illustrated in Fig. |2(b)| The V — A 
bias in the W decay angle causes leptons at high rapidity to 
originate primarily from W bosons produced in the opposite 
hemisphere. 

W bosons at the Tevatron are primarily produced from 
the valence quarks in the proton and the anti-quarks in the 
anti -proton and rarely from sea quarks simply because W pro- 
duction requires at least one moderately high x parton to be 
involved in the collision. At very large forward or backward 
rapidities where one very high x parton must participate in 
the production, the production probability from the sea quarks 
nearly vanishes. Understanding of the sea quark contribu- 
tion is important to exactly know the decay angle distributions 
from the V — A structure because W production by sea anti- 
quarks will result in the opposite W polarization from valence 
quark production. 

We use a Monte Carlo simulation with NLO QCD cor- 
rections to determine the production probability with sea 
quarks by identifying initiating quarks as a function of yw- 
We verify the expected angular distribution of (1 ± cos8*) 2 
from production of W with quarks in the proton and the 
opposite distribution with anti-quarks in the proton. For ex- 
ample, in Fig. |4(a)| we show the cos6>* distributions of e + in 
the W + rest frame for the case when a quark from the proton 
and an anti-quark from the anti-proton form the W + (labeled 
"quark") and the case when an anti-quark from the proton 
and a quark from the anti-proton form the W + (labeled "anti- 
quark"). The ratio of quark (proton) and anti-quark (proton) 
induced W production, therefore, determines the angular de- 
cay distribution. In the simulation, we measure the fraction of 
quark and anti-quark contributions, and parameterize the an- 
gular distributions for yw and the W transverse momentum, 
Pr . We find an empirical functional form that fits the data, 

P±{cos9*,y w ,p¥) = {lTco S e*) 2 +Q{yw,P^){l±cos9*) 2 1 

(6) 

QiVWiPr) = /(^) e -Wp?>^- 2 +o.o5Him- 3 |]. (7) 
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FIG. 4: (a) The cos#* distributions of e + in the W + rest frame, 
averaged over all produced W + . The curve labeled "quark" shows 
the case when a quark from the proton and anti-quark from the anti- 
proton form the W + . The curve labeled "anti-quark" shows the op- 
posite case, when an anti-quark from the proton and a quark from 
the anti-proton form the W + . (b) The dependence of the ratio of 
"anti-quark" (g) and "quark" (q) contributions to the overall W de- 
cay angle distribution, Q(yw,Pr ), as a function of W rapidity and 
px of the W. 

The parameters f(P^) and g(Pjf) are 

!{Pt) = 0.2811£(P^,/i = 21.7GeV,CT = 9.458GeV) 

+0.2185 e (- ao4433GcV " p "\ 
g(P^) = 0.2085 + 0.0074GcV~ 1 P| l/ 

-5.051 x 10 _5 GeV -2 P^ 2 

+ 1.180 x 10" 7 GcV~ 3 P^ 3 , (8) 

where £(x, fx, er) is the Landau distribution with most prob- 
able value fx and the RMS er. The first term of Eq. [6] corre- 
sponds to the contribution from quarks in the proton and the 
second term from anti-quarks in the proton. The parameteri- 
zation, Q(yw ,Pt), the ra tio of the two angular distributions 
as a function of the W rapidity and pj? , is obtained from the 
fit to the distribution in Fig. |4(b)| 

A second relevant factor in the selection among the two W 
rapidity solutions is the W differential cross-section as a func- 
tion of yw, dcr^ /dyw- The W boson production decreases 
sharply beyond \yw\ > 2 because of the scarcity of high x 
quarks. For instance, if one of the two possible solutions falls 
in the central region of rapidity and the other has \yw\ > 2, 
the former should receive more weight as the latter is very 
unlikely to be produced. 

The information used to select among the two solutions can 
be represented by a weighting factor for each rapidity solution 
and charge, wf 2 > can t> e represented as 

± _ P ± (cos9* ]2 ,yi,2,P?' , )T ± (vi,2) 

!> 2 ~~ P± (cosBl ,y i )<r± (yi)+P± (cos8' 2 ,y 2 ,p™ )<r± (y 2 ) ' K ' 

where the ± signs indicate the W boson charge and indices of 
1, 2 are for the two W rapidity solutions. 

In our analysis, we include kinematic cuts for detecting 
charged leptons. For W — > ev event selection, we apply 



|^ ab | < 2.8, E% > 25GeV, and E T > 25GcV. We also 
consider a multiplicative correction factor for the detector ac- 
ceptance and event migration from smearing effects as shown 
in Figs. |5(a)| and |5(b)| In order to study smearing effects, 
we use the fact that the energy resolutions in the electromag- 
netic calorimeter of the Collider Dectector at Fermilab (CDF) 
are (central calorimeter) and 16% /VE © 1% (the 

end plug calorimeter) and in the hadronic calorimeter are 
75%/VE (central) and 80%/VE®5% (the end plug) |J. We 
randomly smear the electron and recoil hadronic energies in 
simulated events with a Gaussian distribution modeling their 
uncertainties prior to making the selection above. The cor- 
rection factors are determined using a Monte Carlo program 
which includes both a model of the process under study as 
well as a simulation of the measuring apparatus. In Eq. [9] 
the weighting factor depends primarily on the W + and W~ 
cross-sections, but does have some weak dependence on the 
assumed W charge asymmetry, and thus the correction factors 
can be biased by computing the factors with different Monte 
Carlo models. Therefore, this method requires us to iterate 
the procedure to eliminate our measurement's dependence on 
the input asymmetry. In order to confirm our analysis tech- 
nique and take into account the bias from physics input vari- 
ables (such as the charge asymmetry itself, the total differen- 
tial cross-section and the angular distribution) we have studied 
the W charge asymmetry measurement with different Monte 
Carlo models and evaluated systematic uncertainties, which 
are described in the next section. 



m. SYSTEMATIC UNCERTAINTIES 

We consider potentially significant sources of systematic 
uncertainty on the W charge asymmetry measurement from 
the assumed parton distributions, the detector resolutions and 
misidentifications and backgrounds. Input PDFs are used to 
determine the parameters of the weighting factor, and may af- 
fect the final result. The detector resolutions affect the W 
rapidity reconstruction due to uncertainties in the calorime- 
ter energy scale and its energy resolution, and the missing 
transverse energy scale also has a significant uncertainty from 
the W boson recoil energy scale. Finally, the detector may 
misidentify the charge, especially from leptons at high \yw\, 
and there are backgrounds to W — ► ev at the Tevatron. 

The uncertainties on the weighting factor (Eq.|9]l arise from 
uncertainties on the momentum distribution of quarks and glu- 
ons in the proton modeled with the PDF sets used. The choice 
of PDF set has an effect on the shape of the dcr ± /dyw distri- 
bution as well as on the ratio of quark and anti-quark in the 
angular decay distribution. We use the CTEQ6 error PDF 
sets lioll and re-determine the da^ /dyw production cross 
section and the angular distribution of (1 ± cosQ*) 2 for each 
error PDF set. We evaluate the uncertainty on the W charge 
asymmetry by checking the deviation of the asymmetry val- 
ues based on each calculation from the central value obtained 
using the best-fitted PDF set. 

We also consider several experimental sources of system- 
atic uncertainty. The scale and resolution of the electromag- 
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FIG. 5: (a) Comparison of the simulated W boson rapidity with 
the generated W rapidity. The simulated rapidity is reconstructed 
by using the weighting factor and the smearing effect is considered, 
(b) Acceptance distribution for W — > ev events as a function of yw- 
The solid (open) circles represent events. 



netic calorimeter energy and the missing transverse energy 
(^t) can change the W rapidity and thus the asymmetry mea- 
surement. We use the energy uncertainties measured in fiUl . 
where the electromagnetic (EM) calorimeter energy scale and 
resolution was tuned in the simulation so as to fit to data. The 
uncertainty on the energy scale and resolution was measured 
to be 0.3% (scale), 1.5% (resolution) for central electron and 
0.6% (scale), 1.1% (resolution) for plug electron. These val- 
ues correspond to a 3 a variation. The asymmetry uncertain- 
ties are estimated as the changes in the measured asymmetry 
when the energy scale and resolution are changed between its 
default and the ±3tr value. The missing transverse energy 
($ T ) in our W — > ev sample is determined by the assumption 
that the vector sum of all transverse energy should be zero. 
Since hadronic transverse energy is due to the W boson re- 
coil energy, we consider the transverse recoil energy, which 
is affected by multiple interactions in the event. The uncer- 
tainty on the transverse recoil energy scale is 2% (3a) frill . 
The charge misidentification rate and background estimates 
are crucial for the charge asymmetry measurement since both 
can directly change the measurement. We estimate these un- 
certainties using the charge fake rates (CFR) and background 
fractions (BKG) from the previous W lepton charge asymme- 
try result from CDF The charge fake rate is about 0.01 
for \rj\ < 1.5 and 0.04 for \r/\ > 2.0. The upper bound on 
the background fraction is 2% for < 1.0 and increases 



TABLE I: Systematic and statistical uncertainties assuming an inte- 
grated luminosity of for the W production charge asymmetry 
in W rapidity bins. 
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to about 15% for |?7| > 2.0. We also investigate sources of 
any charge bias and rj dependence in the kinematic and geo- 
metrical acceptance of the event. An uncorrected acceptance 
shift of 3% central and forward electrons and 5% far forward 
electrons ( | ^7 1 > 2.4) based on measurements of Z — > e + e~ 
data lfl2ll are taken to address the effects of systematic on W 
charge asymmetry measurement. 

Table Q] summarizes the systematic uncertainties on the W 
boson production charge asymmetry for rapidities \yw\ < 
3.0. We compare the expected statistical uncertainty obtained 
by assuming an analysis using an integrated luminosity of 1 
fb _1 , where we also extrapolate the expected statistical uncer- 
tainty from the number of events from the previous W lepton 
charge asymmetry result of CDF with 0.2fb _1 JH. 



IV. RESULTS 

We compare the expected statistical uncertainties in 1.0 
fb _1 of data at the Tevatron with the uncertainties coming 
from parton distribution functions (PDFs) using CTEQ6M in 
Fig-|6] In particular, we notice that at high rapidities {\yw\ > 
1.4) there is a large difference in the precision with which 
the as yet unmeasured W production asymmetry and the pre- 
viously measured asymmetry from the decay leptons scaled 
to 1.0 fb _1 of integrated luminosity are known. The to- 
tal systematic and statistical uncertainties on the W produc- 
tion charge asymmetry measurement is shown in Fig. [7] with 
the uncertainties coming from parton distribution functions 
(PDFs) using CTEQ6M. Since the systematic uncertainty esti- 
mates, as summarized in Table |U are lower than the statistical 
error, a direct measurement of the W charge asymmetry with 
this method should significantly improve parameterizations of 
the PDFs. 

In this paper, we present a study of the W boson produc- 
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tion charge asymmetry with the W decaying leptonically to 
an electron and neutrino in pp collisions at the Tevatron. We 
propose a new analysis technique which resolves the ambigu- 
ity in the neutrino longitudinal momentum, using a realistic 
Monte Carlo simulation. We show that the W charge asym- 
metry can be directly measured at the Tevatron. We conclude 
that by measuring the W production charge asymmetry with 
reconstructed W rapidity, the result should be one of the best 
determinations of the proton d/u momentum ratio, and play 
an important role in global PDF fits. 
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FIG. 6: Comparison of statistical uncertainties expected from this 
analysis in 1.0 fb _1 with those from the uncertainties from CTEQ6M 
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FIG. 7: Total systematic uncertainty estimate and statistical uncer- 
tainty of this method compared with the current uncertainty on the 
W charge asymmetry from the CTEQ6 PDFs. 
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